Stripes and superconducting pairing in the t-J model with Coulomb interactions 
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We study the competition between long- and short-range interactions among charge carriers in 
strongly-correlated electronic systems employing a new method which combines the density-matrix 
renormalization-group technique with a self-consistent treatment of the long-range interactions. 
We apply the method to an extended t-J model which exhibits "stripe" order. The Coulomb 
interactions, while not destroying stripes, induce large transverse stripe fluctuations with associated 
charge delocalization. This leads to a substantial Coulomb-repulsion-induced enhancement of long- 
range superconducting pair-field correlations. 



PACS numbers: 71.10.-w, 74.20.-z, 71.45.Lr 
I. INTRODUCTION 

Much of the interesting physics of the high- 
temperature superconductors, in particular that related 
to the "mechanism" of high temperature superconductiv- 
ity, is moderately local, involving physics on the length 
scale of the superconducting coherence length, ^o- Since 
^0 is typically a few lattice constants, this would seem 
to indicate that numerical solutions of model problems 
on clusters with as few as 50-100 sites should be able to 
provide considerable insight concerning these problems, 
even though results in this range are manifestly sensitive 
to the choice of boundary conditions and other finite-size 
effects. Such studies can also serve as important tests of 
the predictions of analytic theories. 

Studies of t-J systems have, indeed, provided strong 
evidence of a universal and robust d-wave character of 
local pairing correlations^, and of a strong clustering ten- 
dency of holes, which might either lead to "stripe" {i.e. 
unidirectional charge density order) formationcl, or phase 
separations. Both of thesej-features were in fact antic- 
ipated by analytic theoriesclla. However, these studies 
for the t — J and parent Hubbard models have failed to 
find compelling evidence of the strong superconducting 
correlationsH needed to understand high temperature su- 
perconductivity. Moreover, many features of the results, 
especially with regards to stripes, appear very sensitive 
to small chaiiges in the model, e.g. the shape and size 
of the clusteiO, whether or not a small second neighbor 
hopping t' is included or notS, etc: 

All these calculations omit the long-range part of the 
Coulomb interaction because it is difficult to treat using 
any of the standard numerical methods. However, since 
in the high-Tc cuprates, where interactions are gener- 
ally conceded to be strong, the (inter-band) screening 
is semiconductor- and not metallic-like, so there is no a 
priori justification for neglecting the long-range part. In 
addition, in the case of stripes, charge inhomogeneities 
or phase separation, longer-range Coulomb interactions 
are clearly important, a point which has previousljipbeen 
addressed with various mean-field approximationsaEllj 



In this paper, we present a new computational method 
for studying the ground-state properties of electrons with 
strong short- and long-range interactions on fairly large 
finite systems. The method, which may be termed a 
"density-functional DMRG" , uses numerically very ac- 
curate DMRG methods to treat the short-range part of 
the interactions. The long-range piece is taken into ac- 
count within the Hartree approximation, which becomes 
exact in the long-distance limit, and, as we show below 
and in Sec. turns out to work well already for short 
distances. An important point to note is that, through 
the self-consistency requirement, the Hartree potential 
accounts for screening effects in a similar spirit as density- 
functional theory. 

We have studied x 4 t-J ladders and cylinders, with 
hole densities per site rih — 1/9, 1/8, and 1/7, and 
N = 27, 18, 16, and 14. Our principal findings, as sum- 
marized in the figures, are: 1) Charge stripe formation 
is robust to the inclusion of Coulomb interactions of rea- 
sonable strength, although the associated charge density 
modulations have their magnitude somewhat reduced. 
However, "spin stripe" correlations (i.e. spin density 
modulations which suffer a tt phase shift across the charge 
stripe resulting in a spatial period twice that of the charge 
modulations) are prominent on cylinders, but very weak 
on ladders; they are slightly enhanced by the Coulomb 
interactions. 2) The inclusion of Coulomb interactions 
strongly enhances the superconducting pair-field correla- 
tions at the longest distances accessible in these calcula- 
tions. This is our most striking result. 3) We present evi- 
dence that stripe formation does not suppress local super- 
conducting pairing. On the other hand, rigid ftf fpe or- 
dering competes with long-range phase orderingl3^B. The 
enhanced superconducting correlations at long-distances 
produced by Coulomb interactions are, thus, tentatively 
associated with the enhanced pair tunneling between 
stripes produced by the increased stripe fluctuations, in 
agreement with phenomenological argumentaHI. 

This paper is organised as follows. In Sec. || we intro- 
duce our model and the combined DMRG-Hartree tech- 



nique to deal with it. In Sec. HI we present our results 
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for the charge and spin densities as well as for the pairing 
susceptibility. Sec. IV discusses our conclusions. For the 
sake of completeness, a discussion of the validity of the 
Hartree approximation is given in Sec. 



II. MODEL AND TECHNIQUE 



, usins-only a fraction of all the possible states of the 
systemEj. We have typically kept around 800-1000 states 
in the last iterations of the calculation, which results in 
a maximum discarded weight of the order 10~^. We use 
systems with open (and cylindrical) boundary conditions 
chosen not to frustrate the domain walls. 



The Hamiltonian of the t-J model with nearest- 
neighbor hopping t and exchange interaction J plus long- 
range Coulomb interactions, which operates in the sub- 
space of no doubly occupied sites, is given by 

H=-t J2 (cLcr'.+h.C.) (1) 
(r,r')cr 

+ J^(Sr-Sr,-^)+5]l^Coui(rK, 

{r,r') r 

where (r,r') are nearest-neighbor sites, a is the spin in- 
dex, is the electron creation operator, S is the spin 
operator, and Ur — cl^c^a- The long-range part of 
the interactioi£3 



T^coui(r) =Vo J2 



ri(r') 



(2) 



is treated in the self-consistent Hartree approximation, 
whereby the density operator is replaced with its 
ground-state expectation value n(r) =< > and n 
is the uniform positive background charge-density. The 



Coulomb prefactor Vq is given by Vq 



where a is the lattice constant, and r is the coordinate 
of a lattice site in units of a. This long-range poten- 
tial is screened by a background dielectric constant e, 
given both by electronic interband and phonon contri- 
butions, which we take to be e — 8.5, in reasonable 
accordance with , cluster calculations for Coulomb ma- 
trix elementg^nrl and with Quantum-Monte-Carlo simu- 
lationsM. Thus, for a « 3.4 A, the Coulomb prefactor in 
is Vq ~ i, being of the same order of magnitude as 
the kinetic energy. 

Since Vcoui and the density n(r) depend on each other, 
one needs a self-consistent solution of (0) -(H) . Our the- 
oreticaJ treatment is a combination of the DMRG tech- 
niqueEj, which accurately accounts for short-range inter- 
actions, and a Hartree treatment of long-range interac- 
tions. In the spirit of "density-functional theory", we 
iteratively solve (|l]) and (Q) as follows: in the first step 
we set Vcoui — and perform a DMRG calculation. This 
gives rise to a density profile {n{r)} and, via (||) , to a 
new potential Vcoui, which enters the next step DMRG 
calculation as an additional on-site potential. This proee.- 
dure should be repeated iteratively up to convergencellj. 
The fact that the Hartree approximation is appropriate 
for the long-range Coulomb part is discussed in detail in 
Sec. 0. 

Within a given loop, the goal of the DMRG simulation 
is to find iteratively an eigenstate of the Hamiltonian p) 



III. RESULTS 

In the end, we are able to compute ground-state ener- 
gies for given quantum numbers and ground-state cor- 
relation functions. Because we do not have excited- 
state data, we are unable to compute dynamic suscep- 
tibilities. In the figures, we present representative re- 
sults for various ground-state correlation functions. So 
as to make our principal findings clear, we have av- 
eraged these quantities over the transverse coordinate 
of the ladder. Thus, we define the average hole and 
spin density as a function of position along the lad- 
der < a; < iV as p{x) = J2t/=i[^ ~ n{x,y)]/4:, and 
Szix) = E^=i(-1)"+^ < S'ix,y) > /4. Note that 
translational (in the x direction) and spin rotational 
symmetry are explicitly broken in these calculations by 
the ladder ends,^ themselves, and (following White and 
Scalapino (WS)B) by an applied staggerejd-|Zeeman field 
of magnitude h — O.lt on the ladder endsEOn 

To probe superconductivity, we have computed the 
ground-state pair-field correlation function 



D{x) 



N 



At 



N 



(3) 



Heij£_At(a;) creates a (1^2 _y2 -like pair around the (x, 2) 



sitell3. We have explored the dependence of our results 
on parameters to some extent, but in all the figures have 
adopted a conventional value J /t — 0.35. 

Consider first the results of the DMRG calculations 
without the long-range Coulomb potential. Results for 
the spin and charge density are shown as the dashed lines 
in Figs. |l| and || and for the pairing susceptibility D{x) 
by the dashed line in Fig. 

In the low-doping case, [uh = 1/9), depicted in Fig. 
charge stripe order is clearly seen. Stripes form in the 
t-J model so as to satisfy the competing requirements 
of minimizing the kinetic energy of the doped holes and 
minimizing ±lic disturbance of the background exchange 
interactionstlu. In order to distinguish between stripes 
and ordinary Fricdel oscillations, we have compared the 
hole-density profile for different length ladders (18x4 and 
27 X 4) at the same doping. As one can see from Figs. ^ 
and Figs. ||d for the cases with and without Coulomb 
interaction, respectively, the amplitude in the center of 
the system is essentially independent of the system size, 
while Friedel oscillations shojJ^ decay as a function of the 
distance from the boundaryQM. As seen in Fig. |l|b, any tt 
phase shift in the AF order (spin stripes) is quite weak in 
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the present calculation, although it is slightly enhanced 
by the Coulomb interaction. In accordance with Refs.EI 
for the bare t ^ J model, spin stripe order is stronger 
in the case of cylindrical boundary conditions. There 
is roughly one hole per two stripe unit cells, which was 
taken by WSo as evidence that the t — J model favors 
stripes with a minimum energy for a linear charge deiisity 
of A '--^ 0.5. This is in agreement with experimentscj E3, 
which find strifies pSsath A « 0.5 at hole dopings smaller 
thann;,,, «l/8m. 




FIG. 1. (a) Hole density p{x) and (b) staggered spin den- 
sity 5*2(3;) as a function of position x along a 18 x 4[4e,dder 
(i.e., with open boundary conditions in the y direction)ll3 with 
8 holes ijih ~ 1/9). The results with nonvanishing Coulomb 
prefactor (Vb = t, solid lines) are compared with those with- 
out Coulomb interaction (Vb = 0, dashed lines). A compari- 
son of p{x) in a longer system (12 holes in 27 x 4 ) with the 
same = 1/9 (dashed-dotted line in both cases) is shown in 
(c) and (d) respectively with and without Coulomb interac- 
tions. 




FIG. 2. Hole density p{x) for 8 holes in a 16 x 4 (a) and 
a 14 X 4 (b) ladder. Conventions are as in Fig. 0. 

For rifi > 1/8 the density- wave structure seen in Fig. |^ 
is less clear. Were the stripes to retain their integrity at 
these higher hole concentrations, they would be forced 
very closer together, at considerable cost in energy. How- 
ever, the clustering tendency of holes is still apparent in 
the ground-state charge distribution, which is suggestive 
of two hole-rich puddles, each with four holes. 

As discussed by WS, the hole clusters locally share a 
number of features with the two-hole pair state, which 
accounts for the fact that the energy per hole for a do- 
main wall is close to the energy per hole for a paiin. 
This is suggestive of a coniHetition between stripe stabil- 
ity and superconductivityli3. Such a competition has al- 
ready been demonstrated in a model which includes next- 
nearest- neighbor hoppings t' by WS. For large enough 
\t'\, the domain walls "evaporate" into quasiparticles 
{t' < 0) without significant pairing correlations or into 
pairs (t' > Op. It has long been clear that stripe 
formation-suppresses long-range superconducting phase 
coherence^, as is clear from the rapid falloff with dis- 
tance of the pair-field correlator in Fig. ^ (dashed line). 

We now turn to our results in the presence of the 
Coulomb interactions. Unscreened Coulomb interactions 
(e = 1, i. e. Vb — 8.5t) are so strong that they en- 
tirely dominate the physics, destroying all clustering or 
pairing tendencies of the system. However, for the phys- 
ically relevant case, e ~ 8.5, the results are much more 
interesting. 

In the hghtly doped case, Uh = 1/9, shown in Fig. [|, 
the stripe structure is essentially unchanged by Coulomb 
interactions, although the amplitude of the charge modu- 
lations is suppressed (by roughly a factor of 1.5), and the 
anti-phase character of the spin correlations is slightly 
enhanced. We interpret this as meaning that the stripe 
order is robust, but that the Coulomb interactions en- 
hance the transverse stripe fluctuations. 
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FIG. 3. (a) The long-distance part of the pair-field corre- 
lation D{x) for 8 holes in a 18x4 ladder without 1/coui (dashed 
line) and with Vcoui (solid line), (b) The ratio between D{x) 
with Vcoui and without Vcoui for A'^ x 4 ladders with 8 holes 
and iV = 18|-(|iashed) iV = 16 (dotted), and N = 27 with 12 
holes (solid)ta. 

The most dramatic effect of the Coulomb interactions 
is the strong enhancement of the pair-field correlations 
shown in Fig. ^. From Fig.^ (in which the short-distance 
data is off scale) one might conclude that the Coulomb 
interactions primarily increase the overall magnitude of 
the pair correlations. On the other hand, from the ratio 
between the two functions, which is displayed in Fig. ^ 
over the whole range of one can see that it is only the 
long-distance part that is enhanced. In order to sort out 
boundary effects, which are probably responsible for part 
of the strong increase of the ratio at the longest distances 
available, we have also treated a larger (27 x 4) ladder. 
Althouph the calculation of D{x) is less accurate for this 
systemES, we can still draw some conclusions. As one 
can see from the figure, the ratio of D{x) shows oscilla- 
tions with twice the stripe periodicity, whose envelope is 
clearly increasing with distance, even far away from the 
boundaries. 

The dependence of D{x) on distance, seen in Fig. ||a, 
is altered from rapidly decreasing in the absence of 
Coulomb interactions, to much slower distance depen- 
dence in the presence of Coulomb interactions. This re- 
sult supports the idea that while the longer-range phase 
coherence (or, in other words, pair delocalization) is in- 
hibited by rigid stripe order, as obtained in the pure t — J 
model, stripe fluctuations, induced by the Coulomb inter- 



action, permit the pairs to tunnel from stripe to stripe. 
It should be pointed out that this conclusion may not be 
generic and may depend on the stripes stiffness. If stripes 
are intrinsically weak, we expect their fluctuations to be 
strong without the need of Coulomb interaction. In this 
case, the pair-breaking effects of the Coulomb interaction 
would probably obtain the opposite effect and suppress 
pairing correlations. It could also be the case that in our 
relatively narrow systems stripes are particularly stiff, 
due to the fact that they cannot meander effectively. 

To further corroborate this interpretation, we have 
computed the (finite-size) spin gap with and without 
Coulomb interactions for the iV = 14 system (with 
/i = 0). This is As ~ 0.16^ without Coulomb interac- 
tion, and As ^^.I2t with Vq = t. Since we know from 
the work of WSEI, that the local stripe energetics is essen- 
tially dominated by the short-range pair binding (which 
sets a scale for the spin gap as an 5 = 1 excitation of 
the pair) and since wc find only slightly reduced in 
the presence of Coulomb interaction, it is reasonable to 
assume that the local pairing itself is still due to the 
short-range t — J physics. 

The effects of Coulomb interactions on the charge- 
density profile are even stronger (and more complex) in 
the more heavily doped systems as shown for iV = 16 
(uh = 1/8) and iV = 14 (uh = 1/7) in Fig. |. In 
these systems, even the period of the stripe array is al- 
tered by the Coulomb interactions. In particular, the 
two hole puddles seen as minima in the electronic den- 
sity in the absence of Coulomb interactions are broken 
apart into structures that look somewhat more like the 
four-stripe pattern seen in the N = 18 ladder. This result 
is an example of a period selection caused by a Coulomh- 
frustrated tendency to clustering, or phase separationQ. 
In this case, we expect particularly large fluctuations be- 
tween the competing configurations to enhance the delo- 
calization of pairs between stripes. 

IV. CONCLUSIONS 

In summary, in this paper we present a ncijz_approach 
capable of bridging the gap between workBEj arguing 
that stripes are due to a delicate balance between kinetifi 
and long-range Coulomb interaction energy, and workcl 
support the idea that short-range interactions alone can 
lead to stripe formation. Based on a new technique 
( "density- functional" DMRG), we present numerical re- 
sults, supporting the view that — while short-range "t- J- 
like" interactions locally bind holes into pairs — it is the 
long-range Coulomb interaction which induces their de- 
localization accompanied with substantial enhancement 
of superconducting pairing correlations. Moreover, the 
Coulomb-induced stripe fluctuations suppress the mag- 
nitude of charge-density wave order, but can actually 
slightly enhance the ( "anti-phase" ) spin-density wave 
correlations with twice the wavelength. 
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APPENDIX A: VALIDITY OF THE HARTREE 
APPROXIMATION 

The Coulomb interaction part of the Hamiltonian 
reads 

^=^E Ur^-rnr', (Al) 

where is the density operator. This can be rewritten 
as 

y , ^° < nr > nr' (A2) 
^-^ r — r 

+ n / . 1 77 ("r- <nr>) [rir' - < Ur' >) 

2 ^ — ' r — r' 

-|-constant , 

where < nr > is the ground-state expe ctat ion value of rir- 
Here, the first term on the r.h.s. of ( [A^ ) (let's us call 
it Wh) corresponds to the Hartree approximation, while 
the second (Wc) gives its correction. Of course, when the 
expectation value of M^c vanishes, the Hartree approxima- 
tion becomes exact (for ground-state properties). Since, 
by definition, g{r,r') = <^^>''"<^^,> , < Wc > can be also 
written as 

<Wc>^IY. T^^\ <nr><nr<> (5(r,r') - 1) , 

(A3) 

and, thus, the correction to Hartree vanishes when 
g(r,r') = 1. 

In order to measure the accuracy of a mean-field ap- 
proximation for a generic operator 0(r) O(r') in the 
Hamiltonian, one should evaluate its fluctuations 

Ao(r,r') ^ ((6(r)- < 0(r) >) (O(r')- < O(r') >)) . 

(A4) 

Since these terms gives the corrections to the mean-field 
approximation (cf. ( [A3| ) ), a small value for Aq means 
that the mean-field approximation for O is accurate. 



In our calculations, the density fluctuations A„(r, r') 
are typically about 20 — 30 times smaller than the spin 
fluctuations As(r,r'), for nearest-neighbor r and r', 
where both fluctuations are largest (for r ^ r'). This 
is the reason for using the Hartree approximation for the 
Coulomb part and to treat the "J" part exactly. Sim- 
ilarly, on the same site, A„(r, r) « 0.1, which is of the 
same order as As(r, r'), which, in turn, justifies an exact 
treatment of the on-site interaction. 

Notice that the inhomogeneous state brought about 
by the open boundary conditions is decisive in making 
the Hartree approximation work better. Indeed, in a ho- 
mogeneous system (obtained by periodic b.c.) < nr > 
would be constant, and the oscillations in the charge will 
be all shifted to the fluctuations 71^— < > , making 
the correction A„(r, r') large. On the other hand, for 
open boundary conditions, part of the charge oscillations 
are taken care of by the mean values < rir > making 
Ur— < Ur > (and thus A„(r,r')) smaller. 
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